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CONVERSION 

FACTORS,  U.  S.  CUSTOMARY 
UNITS  OF  MEASUREMENT 

TO  METRIC  (SI) 

U.  S.  customary  units 

of  measurement  used  in  this  report  can  be  con- 

verted  to  metric  (SI) 

units  as  follows: 

Multiply 

By 

To  Obtain 

inches 

2.5*4 

centimeters 

feet 

O.30I48 

meters 

square  inches 

6.^516 

square  centimeters 

pounds  (mass) 

0.145359237 

kilograms 

pounds  (mass)  per 
cubic  foot 

16.0185 

kilograms  per  cubic  meter 

slugs  per  cubic  foot 

515.3788 

kilograms  per  cubic  meter 

pounds  (force)  per 
square  inch 

68914.757 

pascals 

pounds  (force)  per 
square  foot 

U7.88026 

pascals 

feet  per  second 

0.30148 

meters  per  second 
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AN  ANALYTICAL  STUDY  OF  PROJECTILE  PENETRATION  INTO  ROCK 


CHAPTER  1 
INTRODUCTION 


1.1  SETTING  THE  STAGE 

Penetration  of  projectiles  into  earth  media  is  a subject  of  great 
interest  due  to  the  increasing  number  of  potential  applications.  With 
respect  to  military  applications,  in  w)  -h  the  projectiles  may  range 
from  small-caliber  munitions  and  shell  fragments  to  bombs  and  missiles 
or  specially  designed  earth  penetrating  weapons  (EPW),^  interest  arises 
from  . oth  offensive  and  defensive  considerations.  Typical  applications 
include: 

1.  Destruction  of  underground  targets  by  penetrating  bombs  or 
missiles  which  explode  near  or  within  the  target. 

2.  Destruction  or  disablement  of  surface  targets  due  to  the  cra- 
tering, ejecta,  or  direct-induced  ground  shock  phenomena  associated 
with  a nuclear  or  high  explosive  detonation  buried  nearby  at  a shallow- 
to-optimum  depth. 

3.  Implantation  of  mines  and  sensors. 

1*.  Design  and  construction  of  undei  ground  protective  structures 
to  resist  penetrating  weapons  effects. 

McNeill  (Reference  i)  describes  the  potential  use  of  instrumented  pene- 
trators  for  rapid  mapping  of  soil  and  rock  profiles.  Typical  situations 
might  be  the  investigation  of  alternate  construction  sites  and  alternate 
routes  for  trenches,  canals,  and  tunnels  and  the  study  of  inaccessible 
and  remote  areas. 

This  report  presents  the  results  of  a parametric  study  of  pro- 
jectile penetration  into  rock.  Rock  penetration  has  not  been  investi- 
gated in  laboratory  and  in  situ  test  programs  as  extensively  as  has 


For  convenience,  symbols  and  unusual  abbreviations  are  listed  and 
defined  in  the  Notation  (Appendix  A). 
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penetration  into  soils.  As  military  penetration  interest  extends  to  in- 
creasingly greater  depths  of  penetration,  the  probability  of  encounter- 
ing rock  increases  rapidly;  also,  there  is  a significant  percentage  of 
the  earth's  surface  with  exposed  rock  or  rock  covered  by  only  a thin 
mantle  of  soil.  Thus  penetration  into  rocks  must  at  least  be  considered 
in  the  "trade-off"  studies  involved  in  many  penetrating  weapon  system 
designs. 

In  order  to  effect  penetration  into  rock,  design  structurally 
sound  penetrators  for  use  in  rock,  and  predict  depth  of  penetration,  it 
is  important  that  rock  penetration  be  more  thoroughly  investigated. 
Reference  2 describes  a continuing  program,  initiated  in  1966  by  Sandia 
Laboratories,  on  the  penetration  of  in  situ  rock;  it  is  concluded  in 
this  reference  that  it  is  possible  to  design,  construct,  and  instrument 
projectiles  which  can  survive  the  environment  of  impact  and  penetration 
into  low-  and  medium-strength  rock.  A brief  history  of  rock  penetra- 
tion studies  is  given  in  Reference  2,  and  Reference  3 is  a comprehensive 
state-of-the-art  report  on  the  subject.  From  the  literature  it  is  evi- 
dent that  the  number  of  carefully  documented,  full-scale,  in  situ  pene- 
tration tests  is  very  small.  Thus  th.-  empirical  data  base  available  is 
not  yet  sufficient  to  provide  an  adequate  empirical  basis  for  determin- 
ing the  effects  of  independent  variables  on  the  rock  penetration  process. 
If  a theory  could  be  found  which  would  even  approximately  replicate 
the  limited  data  on  rock  penetration,  it  could  be  used  for  parameter 
studies  of  the  effects  of  its  independent  variables,  and  the  results  of 
this  study  would  then  be  relevant  to  the  performance  of  penetrating 
weapons  impacting  rock.  By  conducting  a systematic  parameter  study 
using  an  analytical  penetration  model  which  does  not  rely  on  empirical 

"constants"  and  which  has  been  successfully  used  in  predicting  depth 

2 

of  penetration  for  a number  of  shallow  penetration  events  into  a variety 
of  materials,  including  rock  and  concrete,  it  is  hoped  that  this  report 
will  allow  the  reader  to  gain  a qualitative  understanding  of  the  factors 
of  importance  in  rock  penetration,  a feel  for  the  depths  of  penetration 


See  Chapter  U. 
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which  may  be  attainable,  and  an  appreciation  of  the  complexities  and 
uncertainties  involved. 

1.2  PURPOSE 

The  purposes  of  this  report  are  (l)  to  demonstrate  that  a rigid 
body  penetration  theory  developed  from  the  cavity  expansion  theory  (Ref- 
erences I4  through  6)  can  approximately  predict  the  results  of  shallow 
penetration  tests  in  hard  media  and  (2)  to  parametrical] y investigate 
the  effects  of  variations  in  rock  properties,  nose  shape,  sectional  pres- 
sure, and  impact  velocity  on  penetration  depth  and  deceleration-time 
histories  over  a range  which  is  of  potential  interest  in  connection  with 
the  design  of  EPW’s. 

1.3  PENETRATION  PARAMETERS 

Penetration  of  a projectile  into  rock  (or  actually  any  material) 
is  a function  of  parameters  describing  the  impact  geometry,  projectile 
shape  and  material  properties,  and  rock  material  properties.  The  fol- 
lowing paragraphs  describe  the  parameters  involved  in  the  general  pene- 
tration problem  and  then  discuss  the  assumptions  and  specializations 
used  in  this  report. 

1.3.1  Impact  Geometry.  The  geometry  for  a general  impact,  situa- 
tion is  shown  in  Figure  1.1a  and  involves  specification  of  an  angle  of 
obliquity  and  an  angle  of  attack  (also  called  angle  of  yaw);  note  that 
the  projectile  axis  of  symmetry  does  not  necessarily  lie  in  the  plane 
of  incidence,  and  the  projectile  may  or  may  not  be  spinning  about  its 
axis  at  impact.  Treatment  of  this  general  impact  geometry  is  beyond 
the  state  of  the  art,  and  this  report  considers  the  case  of  normal  im- 
pact of  a nonspinning  projectile  (Figure  1.1b)  with  the  trajectory  co- 
incident with  the  projectile  axis  of  symmetry.  Figure  1.1c  defines  the 
depth  of  penetration  as  used  in  this  report. 

1.3.2  Projectile  Shape  and  Material  Properties.  In  general, 
projectiles  of  interest  may  be  considered  to  be  cylindrical,  having  a 
body  length  £ , diameter  D , and  nose  section  which  may  be  blunt 
(right- circular  cylinder),  hemispherical,  conical,  ogival,  or  some 
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combination  of  these  shapes.  The  projectiles  considered  in  this  report 
have  ogival  nose  shapes,  with  the  nose  length  and  '•sharpness”  determined 
by  the  Caliber  Radius  Head  (CRH)  and  the  body  diameter  D as  illus- 
trated in  Figure  1.2.  For  the  normal  penetration  case  considered  in 
this  report,  the  projectile  body  length  t is  not  an  explicit  parameter. 
Also,  the  projectiles  are  considered  to  be  rigid  bodies;  hence  internal 
structure  and  material  properties  are  not  explicitly  considered  (body 
length,  internal  structure,  and  material  density  enter  as  implicit  pa- 
rameters in  the  determination  of  projectile  mass  m). 

1.1*  SCOPE 

The  penetration  model  used  for  the  work  reported  herein  is  pre- 
sented in  Chapter  2.  A discussion  of  a general  rock  classification 
scheme,  which  will  later  be  shown  to  be  useful  for  predicting  penetra- 
tion into  rock,  is  presented  in  Chapter  3,  and  the  parameters  necessary 
for  use  in  the  penetration  model  are  selected.  In  Chapter  1*,  studies 
performed  to  demonstrate  that  the  penetration  model  can  qualitatively 
simulate  many  aspects  of  actual  penetration  tests  inx>  rock  are  pre- 
sented. In  Chapter  5,  a systematic  parameter  study  of  the  independent 
variables  over  the  range  of  interest  in  rock  penetration  is  presented. 
Depths  of  penetration  are  computed  for  a range  of  impact  velocities  for 
projectiles  with  various  weights,  diameters,  and  CRH  values. 
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a.  Geometry  of  generalized  impact  problem. 


b.  Normal  impact.  e.  Definition  of  depth 

of  penetration. 


Figure  1.1  Projectile  impact  geometry. 
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CHAPTER  2 


PENETRATION  MODEL 

2.1  FORCE  LAWS  AND  EMPIRICAL  RELATIONS 

References  1 and  I4  present  histories  of  the  attempts  to  develop 
equations  for  predicting  the  depth  of  projectile  penetration  into  earth 
materials.  The  prevalent  procedures  have  been  to  develop  purely  em- 
pirical equations  relating  the  measured  depths  of  penetration  for  a set 
of  experimental  tests  to  various  test  parameters  or  to  assume  that  the 
penetration  event  is  governed  by  a particular  form  of  force  law.  Using 
initial  and  end  conditions,  the  equation  of  motion  (based  on  the  assumed 
force  law)  may  be  integrated  to  yield  time  histories  and  an  equation 
relating  final  depth  of  penetration  to  initial  velocity.  The  resulting 
equation  will  contain  one  or  more  constants  which  must  be  determined 
from  actual  penetration  data  by  curve-fitting  techniques. 

The  basic  deficiency  of  both  the  purely  empirical  approach  and  the 
force  law  method  is  that  they  both  contain  coefficients  not  uniquely 
defined  in  terms  of  the  physical  properties  of  the  target  materials  and 
the  geometry  of  the  projectile.  Also,  since  the  coefficients  must  be 
evaluated  from  penetration  test  results  for  a given  velocity  range, 
projectile  geometry,  and  target  type,  the  penetration  equations  cannot 
be  used  reliably  outside  the  range  of  test  conditions  used  in  fitting 
the  equations. 

2.2  THEORETICAL  PENETRATION  MODEL 

2.2.1  Original  Formulation.  A theory  of  penetration  of  rigid 
projectiles  into  incompressible  elastic-plastic  targets  was  developed 
by  Goodier  (Reference  5).  The  theory  was  extended  by  Ross  and  Hanagud 
(Reference  6)  to  include  the  compressibility  of  the  target  and  has  been 
used  successfully  in  studies  of  penetration  of  ice,  frozen  soil,  sand, 
and  clay  (References  6,  and  7).  The  theory  has  also  been  used  to  pre- 
dict mine  penetration  into  a variety  of  earth  targets  (Reference  8)  and 
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bomb  penetration  into  concrete  (Reference  9)*^  Treatment  of  penetration 
into  a multilayered  medium  is  presented  in  Reference  10. 

The  penetration  theory  is  based  on  the  theory  of  dynamic  cavity 
expansion  in  locking  elastic-plastic  materials  whose  behavior  can  be 
idealized  by  the  compressibility  and  shear  behavior  illustrated  in  Fig- 
ure 2,1,  in  which  e^  is  the  volumetric  strain  related  to  the  Jump 

from  zero  to  finite  stress  (elastic  region),  e is  the  volumetric 

P 

strain  related  to  the  Jump  from  the  elastic  to  the  plastic  region  of 
the  stress-strain  curve,  E is  the  strain-hardening  modulus  correspond- 

w 

ing  to  the  locked  plastic  region,  0 is  the  strain-hardening  shear  modu- 

v 

lus,  E is  Young's  modulus  of  elasticity  corresponding  to  the  locked 

elastic  region,  and  G is  the  elastic  shear  modulus.  The  resulting 

penetration  equation  is  the  solution  of  the  equation  of  motion  of  a 

rigid  body  subjected  to  a distributed  normal  stress  on  its  nose.  The 

local  normal  stress  is  determined  from  the  theory  of  dynamic  cavity 

expansion  for  an  assumed  velocity  field  in  the  media  adjacent  to  the 

projectile.  The  normal  stress  or  cavity  pressure  is  obtained  from  the 

solution  of  a boundary-value  problem  and  is  expressed  in  terms  of 

parameters  defined  explicitly  in  terns  of  measurable  material  properties 

of  the  target  and  the  mass  and  geometry  of  the  projectile. 

2 

2.2.2  WES  Extensions  of  the  Penetration  Model.  In  the  original 
formulation,  Ross  and  Hanagud  assumed  a hemispherical-nosed  projectile 
and  normal  impact.  WES  is  currently  involved  in  theoretical  efforts  to 
extend  the  penetration  model  to  apply  to  oblique  impact  conditions  as 
well  as  to  conical  and  ogival  nose  geometries  and  layered  targets. 

Since  this  report  considers  ogival  projectiles,  the  penetration  equation 
as  modified  to  apply  to  ogives  is  used.  The  relevant  equations  are 
listed  without  derivation,  since  the  derivations  and  documentation  of 


In  this  reference,  results  of  application  of  this  theory  are  shown 
to  be  in  very  good  agreement  with  empirically  based  concrete  pene- 
tration formulas  within  the  range  of  the  parameters  for  which  the 
formulas  are  considered  accurate. 

The  expression  "penetration  model"  is  used  in  the  reminder  of  this 
report  to  refer  to  the  Ross-Hanagud  penetration  theory  and  any  WES 
modifications  and  extensions  to  the  basic  theory. 
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the  model  will  be  the  subject  of  a future  WES  technical  report. 

The  equation  of  motion  for  a rigid  projectile  of  mass  m and  di- 
ameter D penetrating  in  the  z-direction  (origin  at  ground  surface, 
positive  downward)  is  given  by  the  following: 


m 


£z 

dt2 


'L1  2 at2 


(2.1) 


The  parameters  in  Equation  2.1  are  given  in  terms  of  the  target  material 
properties  by  the  following: 


B1  = 1 - 61/3  (2.2) 

B2  " I " (1  + “p)  61/3  + I fiU/3  (2*3) 

pg  -|E  [1  - exp(-3B)]  - | Y £n  6 + |y  7t2Et  - | Etn  (2.M 

P 

6=1-—  exp(-36)  (2.5) 

PP 

a = 1 - — (2.6) 

P P„ 


n =■  Z h (2-7) 

n=l  n 

(2-8> 

Pp  * PQexp(ep)  (2.9) 


where 

P0  = initial  density  of  target  material,  slugs/cu  ft3 
Y = yield  strength  of  target  material  (shown  in  Figure  2.1),  psf 
E = Young's  modulus  of  elasticity,  corresponding  to  locked  elastic 


region  in  Figure  2.1,  psf  * 3G  (G  * elastic  sheur  modulus) 

E * strain-hardening  modulus,  corresponding  to  locked  plastic 

region  in  Figure  2.1,  psf  = 3G  (G  = strain-hardening  shear 
modulus ) 

* volumetric  strain  related  to  t.h*»  jump  from  zero  to  finite 
stress  (elaslic  region)  in  the  stress-strain  curves  in 
Figure  2.1 

e = volumetric  strain  related  to  the  Jump  from  the  elastic  to 
^ the  plastic  region  of  the  uniaxial  stress  (i.e.,  unconfined 
compression)  stress-strain  curve  in  Figure  2.1b 

The  function  f(t)  is  the  ogive  nose  factor:*4 


f(e) 


- *>c)r  ■> 

■ ? / h 


# + 


3 2 

sin  cos  p 

1 - (l  - e)  cos  ip 


] 


(cos  ip  - 1 + e)  dtj/ 


(2.10) 


where 


1 

2 CRH 


cos  p = 1 - c , or 
4>U)  = tan"1  (-f”") 


(2.11) 


(2.12) 


CRH  is  defined  in  Figure  1.2  and  $ is  the  limiting  cone  half-angle. 

A plot  of  f(e)  versus  CRH  is  given  in  Figure  2.2. 

The  target  material  is  defined  by  six  material  parameters  (e^  , e 

or  p , p , E , E.  , and  Y)  and  the  projectile  by  three  parameters 
p o t< 

(m  , CRH  , and  D).  The  parameters  e.  and  e can  be  obtained  from 
a hydrostatic  compression  test  on  undisturbed  specimens  of  target  ma- 
terial conducted  at  strain  rates  and  carried  to  peak  stress  levels  gen- 
erally anticipated  during  the  penetration  event.  The  parameters  E , 


3 

A table  of  factors  for  converting  U.  S.  customary  units  of  measure- 
^ ment  to  metric  (Si)  units  is  presented  on  page  5. 

R.  S.  Bernard,  private  communication,  January  197^.  Equation  2.10 
represents  a first  attempt  at  incorporating  nose  shape  effects  into 
the  model. 
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, and  Y can  be  determined  from  an  undrained  dynamic  uniaxial  stress 
test  (unconfined  compression)  or  from  a dynamic  triaxial  test  at  a 
confining  pressure  simulating  the  in  situ  stress  state. 

The  constitutive  model  assumed  in  the  cavity  expansion  theory  is 
a relatively  simple  one.  Strength  is  independent  of  stress  level, 
volume  change  can  occur  only  at  the  initiation  of  stress  or  at  the  onset 
of  plastic  behavior,  and  dilatation  in  shear  is  not  possible  in  the 
model.  Such  simplifications  are  at  variance  with  the  observed  behavior 
of  real  soil  and  rock.  Thus,  it  should  be  noted  that  the  material 
properties  used  provide  only  a crude  description  of  the  3tress-strain 
behavior  of  the  material;  and  as  with  any  analytical  technique,  it  is 
within  the  realm  of  possibility  to  obtain  "good  results"  from  a fortui- 
tous selection  of  parameters,  even  though  individually  the  parameters 
might  not  be  the  correct  choice  for  the  particular  event. ^ Their  use 
(as  well  as  the  use  of  the  theory  itself)  can  only  be  Justified  on  the 
basis  of  agreement  of  the  analytical  results  with  available  penetration 
data.  Such  agreement  will  be  demonstrated  later  in  this  report  for  a 
number  of  cases  within  the  range  of  variables  of  relevance  to  EFW  rock 
penetration. 

Equation  2.1  is  easily  amenable  to  an  incremental  solution  in 
terms  of  velocity  by  using  the  initial  and  end  conditions. 


4f-  ■ V , when  z * 0 
dt  o 

dz  „ 

T7-  * 0 , when  z * P 


(2.13) 


where  VQ  is  the  impact  velocity  and  P Is  the  final  depth  of  pene- 
tration. After  n time  steps  At  , the  velocity  V is  given  by 


The  penetration  event,  for  example,  may  produce  a distinctive  stress 
field  in  the  target  which  in  itself  would  establish  the  appropriate 
confining  pressure  for  laboratory  property  determination  (rather  than 
the  natural  in  situ  stress  field). 
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(2.  ill) 


V(nAt) 


where  AV.  is  the  incremental  velocity  change  in  the  k*1*1  time  step 
as  computed  from  Equation  2.1.  The  acceleration  during  the  k^  time 
step  is  AV^/At  , and  the  distance  penetrated  during  the  k time  step 
is  computed  as 


Az  V(kAt)  + V[(k  - l)At]  At  (2.15) 

k 2 

A number  of  small  computer  codes  have  been  used  to  solve  Equa- 
tion 2.1  in  the  manner  described.  Although  it  can  be  solved  in  closed 
form  to  give  a final  depth  of  penetration  in  a homogeneous  half  space, 
solving  the  equation  incrementally  has  a number  of  advantages.  During 
the  early  phase  of  penetration,  until  the  ogival  nose  is  completely 
embedded,  the  effective  cross-sectional  area  of  the  projectile  varies; 
and  by  letting  the  diameter  vary  with  penetration  in  Equation  2.1,  as 
required  to  reflect  the  effective  cross-sectional  area,  the  embedding 
process  is  included  in  the  model.  Also  the  incremental  solution  allows 
the  material  properties  of  the  target  to  be  changed  at  any  time  step, 
which  in  turn  allows  layered  targets  to  be  considered.  The  computer 
codes  may  be  accessed  from  time  sharing  terminals  and  are  easy  and 
inexpensive  to  operate.  On-line  plotters  can  be  used  to  produce  pa- 
rameter plots. 

One  additional  WES  modification  to  the  theory  was  to  allow  p to 

P 

be  changed  with  time  in  the  incremental  method  as  the  cavity  pressure 
drops,  i.e.,  the  complete  target  pressure-volumetric  strain  curve  may 
be  used  instead  of  the  locking  approximation. 
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PRESSURE 


• ELASTIC  RANGE 
O PLASTIC  RANGE 


VOLUMETRIC  STRAIN 


a.  Behavior  in  hydrostatic  b.  Behavior  in  unconfined 

compression.  compression. 

Figure  2.1  Idealized  stress-strain  curves  for  a locking 
elastic-plastic  material. 


Figure  2.2  Ogive  nose  factor. 
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CHAPTER  3 


SELECTION  OF  TARGET  PARAMETERS 

3.1  INTACT  AND  IN  SITU  ROCK  PROPERTIES:  AN  OVERVIEW 

Rock  is  properly  described  as  a nonlinear,  hysteretic  material  and 
in  its  in  situ  state  is  nonhomogeneous , anisotropic,  and  discontinuous. 
Even  the  most  sophisticated  finite  difference  or  finite  element  computer 
code  cannot  treat  the  intrinsic  and  discontinuum  properties  of  the  in 
situ  rock  mass  in  all  their  complexity.  Rock  mechanics,  being  a rela- 
tively new  formal  discipline,  does  not  benefit  from  a large  repertoire 
of  carefully  documented  cases  of  engineering  successes  and  failures. 

Also,  unlike  soil  mechanics,  rock  mechanics  has  not  enjoyed  much  success 
in  using  laboratory  test  results  to  design  and  predict  the  behavior  of 
structures  in  and  on  rock.  In  a fundamental  sense,  the  reasons  for  the 
problems  encountered  in  rock  mechanics  might  be  reduced  to  a problem  in 
probability  and  statistics  (Reference  11 J,1  although  such  an  approach 
may  prove  to  be  of  limited  utility,  since  it  is  difficult  to  conceive 
how  probability  density  functions,  for  example,  could  be  constructed  for 
general  cases. 

Rock  contains  not  only  small-scale  defects  such  as  cracks,  flaws, 
inclusions,  grain  boundaries,  etc.,  but  also  large-scale  inhomogeneities, 
such  as  layers.  Joints,  faults,  fractures,  folds,  etc.,  with  scales 
ranging  from  centimeters  to  hundreds  of  meters.  Also  the  rock  surfaces 
involved  in  these  geologic  structures  may  be  weathered  or  may  be  char- 
acterized as  being  completely  random  with  regard  to  asperities  or  as 
having  undulating  asperities  with  characteristic  "wave  lengths"  and 
"amplitudes"  ranging  from  centimeters  to  hundreds  of  meters.  Thus  the 
reasons  for  the  existence  of  a "size  effect"  problem  in  the  determina- 
tion of  constitutive  properties  become  apparent.  Not  only  are  there 
variations  in  strength  and  deformability  of  intact  rock  specimens  in 


This  reference  contains  a discussion  of  reducing  rock  mechanics  prob- 
lems to  a problem  in  probability  and  statistics  and  also  a comprehen- 
sive bibliography. 
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the  laboratory  as  size  increases  (References  11  P" ' 12),  but  also  the 
mode  of  failure  may  change  (brittle  to  ductile,  for  example).  The  "size 
effect"  problem  becomes  even  more  significant  when  it  is  desired  to  ex- 
trapolate laboratory  data  for  properties  in  situ.  In  situ  tests  have 
been  a tremendous  aid  to  engineers,  but  the  obvious  problem  still  arises, 
as  emphasized  by  Heuze  in  Reference  13,  that  the  size  of  the  in  situ 
test  should  be  representative  of  the  size  of  the  in  situ  prototype. 

Since  in  situ  rock  mass  strength  and  deformability  depend  on  both 
the  intact  properties  and  the  discontinuum  properties,  considerable  ef- 
fort has  been  expended  in  recent  years  to  understand  the  mechanics  of 
Jointed  rock  (References  13  through  15,  for  example).  Also,  progress  has 
been  extensive  in  developing  engineering  oriented  classifications  and 
index  properties  for  intact  and  in  situ  rock  (References  16  through  18). 
The  result  of  these  efforts  is  the  recognition  of  correlations  between 
certain  laboratory  test  results,  between  certain  in  situ  test  results, 
between  laboratory  and  in  situ  test  results,  and  importantly  between 
discontinuum  index  descriptors  and  laboratory  and  in  situ  test  results. 

Obviously,  the  penetration  model  is  intended  to  apply  to  a homo- 
geneous isotropic  continuum.  Thus,  if  it  is  to  be  used  successfully 
in  rock  penetration,  a methodology  must  be  developed  for  selecting 
equivalent  continuum  properties  for  the  rock  target  such  that  the  theory 
gives  answers  which  correlate  with  real  tests  on  the  discontinuum.  A 
methodology  for  this  purpose  is  described  in  this  chapter. 

3.2  USE  OF  IN  SITU  AND  INTACT  ROCK  CLASSIFICATION  SCHEMES 

TO  ESTIMATE  RANGES  OF  CONSTITUTIVE  PROPERTIES 

Reference  17  presents  a classification  scheme  for  intact  rock 
based  on  two  properties  of  engineering  significance  (strength  and  de- 
formability).  Table  3.1,  which  summarizes  the  classification  scheme 
developed  by  Deere  and  Miller  (Reference  17),  evolved  from  numerous 
laboratory  tests  on  many  different  rock  types.  Use  of  the  modulus  ratio 
concept  results  from  an  observed  correlation  between  uniaxial  com- 
pressive strength  and  tangent  modulus  of  elasticity;  Figure  3.1  (from 
Reference  17)  illustrates  the  classification  concept.  Rocks  can  be 
classified  under  this  scheme  from  a knowledge  of  the  values  of  strength 
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and  modulus  determined  in  laboratory  tests.  As  illustrated  in  Fig- 
ure 3.1  for  the  particular  case  of  ignec’is  rocks,  rocks  of  a given  type 
tend  to  cluster  in  a restricted  region  of  the  space.  The  limiting 
values  for  the  average  modulus  ratio  band  were  chosen  to  encompass  the 
majority  of  the  data  for  rocks  with  interlocking  textures.  Thus  some 
of  the  power  and  utility  of  the  classification  scheme  can  be  realized. 

In  many  cases,  knowledge  of  rock  type  and  texture  may  enable  an  "edu- 
cated estimate"  of  engineering  properties  to  be  made;  and  if  rock  type 
and  either  strength  ££  modulus  are  known,  then  the  other  property  can 
be  predicted  with  a much  higher  degree  of  confidence. 

Other  important  correlations  exist  between  the  engineering  proper- 
ties and  various  quantitative  rock  index  properties,  such  as  Schmidt 
hardness.  Shore  hardness,  and  dry  unit  weight.  The  importance  of  these 
correlations  is  that  the  index  tests  can  be  performed  easily,  quickly, 
and  economically  on  rock  cores  in  the  field.  Nomograms  based  on  the 
observed  correlations,  such  as  Figures  3.2  and  3.3,  can  then  be  used  to 
estimate  strength  and  modulus.  Thus  rocks  can  be  classified  according 
to  Table  3.1  based  solely  on  simple  index  tests. 

Reference  18  presents  a systematic  classification  scheme  for  the 
quality  of  in  situ  rock  (i.e.,  the  degree  to  which  the  rock  represents 
a continuum).  The  in  situ  classification  scheme  is  based  on  a quanti- 
tative core  logging  procedure  proposed  by  Deere  in  196U.  The  in  situ 
nature  of  the  rock  mass  is  described  by  its  Rock  Quality  Designation 
(RQD),  which  is  obtained  by  measuring  the  total  length  of  all  un- 
weathered pieces  of  core  greater  than  or  equal  to  U inches  and  then 
dividing  by  the  total  length  of  core  run.  RQD  has  been  shown  to  cor- 
relate with  other  qualitative  descriptions  of  rock  quality  such  as  frac- 

2 

ture  frequency  and  velocity  ratio  (V^/V^)  , where  is  the  in  situ 

and  the  laboratory  (intact)  compressional  wave  velocity.  The  RQD 
is  expressed  as  a percentage,  with  an  RQD  of  100  percent  corresponding 
to  a fracture  frequency  of  0 and  a velocity  index  of  1.  Table  3.2 
summarizes  the  classification  system. 

On  the  basis  of  extensive  in  situ  testing  to  determine  apparent 
elastic  moduli,  such  as  plate  Jack  tests,  radial  Jack  tests,  borehole 
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deformation  tests,  etc.,  various  important  correlations  between  RQD 
and  the  ratio  of  intact  to  in  situ  constitutive  properties  have  been  de- 
termined {References  18  and  19).  Correlations  have  been  established  be- 
tween the  ratio  of  apparent  Young's  moduli  (E^ne^u)  and  the  intact 
specimen  Young's  modulus  (Eintact)  611:1(1  RQD  • Figure  3.1*  shows  one  of 
the  correlations  as  given  in  Reference  20.  The  general  trend  of 
Einsiu/Eintact  versus  RQD  has  ^een  confirmed  by  other  investigators 
in  field  tests,  model  studies,  and  finite  element  studies,  although  a 
straight  line  correlation  should  not  be  expected  and  it  is  probable  that 
the  ratio  approaches  a value  of  0.1-0.15  at  low  values  of  RQD  . This 
suggests  that  an  estimate  of  in  situ  modulus  of  elasticity  can  be  easily 
obtained  from  a knowledge  of  RQD  , Schmidt  hardness,  and  dry  unit  weight 
by  using  correlations  such  as  those  shown  in  Figures  3.3  and  3.1*.  If 
the  intact  modulus  is  known,  only  an  estimate  of  RQD  is  needed  in 
order  to  estimate  the  in  situ  modulus. 

In  many  instances  rock  masses  nay  be  described  as  being  laminated 
or  tri-orthogonally  jointed.  For  these  conditions,  Duncan  and  Goodman 
(Reference  21)  derive  an  expression  relating  in  situ  modulus  of  deforma- 
tion Einsitu  to  the  intact  modulus  Eintact  » normal  Joint  stiffness 

K , and  fracture  spacing  s : 
n 


Jinsitu 


(3.1) 


Jintact 


K s 
n 


By  using  a correlation  by  Deere  (Reference  20)  between  RQD  and  fracture 
frequency  f , where 


(3.2) 


K in  terms  of  E., 


if  it  is  assumed 


it  is  possible  to  determine 

n intact 

that  the  ratio  Einsitu^Eintact  approaches  0.1-0.15  at  RQD  * 0 as 

suggested  previously.  Figure  3.5  shows  the  relation  for  the  ratio 
values  0.1  said  0.15  at  RQD  = 0 . 

The  discussion  thus  far  has  not  referred  to  in  situ  strength,  and 
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indeed  "handy"  correlations  such  as  those  shown  in  Figure  3.5  for  the 
relation  of  in  situ  and  intact  unconfined  compressive  strengths  do  not 
exist.  Thus  it  is  in  general  not  possible  to  deduce  in  situ  strength 
from  a knowledge  of  RQD  and  intact  strength,  for  example,  as  it  was 
for  the  modulus  values  (Figure  3.5).  Much  of  the  problem  arises  from 
misconceptions  regarding  the  nature  of  rock  strength;  as  stated  by 
Hudson  (Reference  22),  properties  such  as  tensile  and  compressive 
strengths  are  experimental,  not  intrinsic,  properties.  Thus  laboratory 
strength  values  depend  on  the  method  of  testing  and  size  of  the  test 
specimen.  Despite  these  concepts,  properties  such  as  tensile  and  com- 
pressive strengths,  yield  strength,  ultimate  strength,  crushing  strength, 
etc.,  are  still  used  and  needed  for  many  types  of  engineering  analyses. 
v Therefore,  caution  is  i*equired  in  specifying  strength  values  and  in 
describing  the  test  method.  Speaking  in  terms  of  a statistically  sig- 
nificant number  of  tests  for  a given  type  of  strength  test,  several 
"values  of  strength"  may  be  observed  as  the  size  of  the  test  specimen 
is  increased,  because  in  some  cases  failure  occurs  due  to  flaws  or 
discontinuities  (References  11,  12,  and  22)  and  the  scale  dimensions  of 
the  various  types  of  discontinuities  vary  by  several  orders  of  magnitude 
(grain  boundaries,  cracks,  layer  boundaries,  joint  block  boundaries, 
faults,  etc.).  Thus  it  is  i nportant  that  strength  data  or  estimates 
apply  to  the  size  scale  of  interest. 

In  order  to  discuss  in  situ  strength  for  projectile  penetration 
calculations,  it  is  necessary  then  to  specify  both  the  type  of  experi- 
mental strength  value  and  the  scale  of  the  test.  For  the  depths  of 
interest  in  projectile  penetration  (<200  feet),  confining  pressures  will 
probably  not  exceed  150  psi.  Thus  triaxial  tests  with  very  low  con- 
fining pressures  or  unconfined  tests  will  suffice.  The  maximum  projec- 
tile body  diameter  considered  in  this  report  is  10  inches,  which  can 
be  considered  as  a representative  size  scale  for  a large  class  of  earth 
penetrators  (Reference  23).  Nearly  all  earth  penetrators  have  a sharply 
pointed  nose  section  (the  projectiles  considered  in  this  report  are  all 
ogives).  Compressive  failure  of  the  rock  occurs  in  the  region  surround- 
ing the  tip  of  the  nose,  and  failure  is  envisioned  to  progress  radially 
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outward  from  the  initial  region  of  failure  as  the  nose  penetrates;  thus 
the  size  scale  of  interest  for  strength  value  considerations  is  of  the 
order  of  or  smaller  than  the  body  diameter.  Based  on  these  considera- 
tions, it  is  postulated  that  strength  values  determined  from  typical 
laboratory-size  test  specimens  are  adequate  for  penetration  calculations 
and  that  corrections  for  the  differences  between  the  appropriate  scale 
in  situ  test  and  the  laboratory  strength  test  are  not  necessary. 

3.3  BASIC  CRITERIA  FOR  THE  SELECTION  OF 

ROCK  TARGET  CONSTITUTIVE  PROPERTIES 

In  Section  2.2.2,  a penetration  model  requiring  six  constitutive 
parameters  for  the  target  is  presented.  The  appropriate  in  situ  Young's 
modulus  E can  be  estimated  from  a value  for  the  intact  (laboratory) 
modulus  and  the  in  situ  rock  quality  RQD  by  the  methodology  presented 
in  Section  3.2.  It  is  fortunate,  for  calculational  purposes,  that  most 
rocks  at  low  confining  pressures  (<1000  psi)  exhibit  brittle  failure  or 
at  least  little  or  no  strain  hardening  (References  2k  and  25).  For 
these  rocks,  the  yield  strength,  ultimate  strength,  compressive  strength, 
and  crushing  strength  can  all  be  considered  the  same  in  unconfined  or 
triaxial  compression  tests  at  low  confining  pressure.  Thus  the  strain- 
hardening Young's  modulus  E^  can  be  set  equal  to  zero,  and  the  value 
for  the  strength  Y can  be  determined  from  unconfined  compression  tests 
on  laboratory-size  specimens.  For  finely  Jointed  and  weathered  rock 
(low  RQD),  it  may  be  necessary  to  apply  reduction  factors  to  the  in- 
tact (laboratory)  value  for  Y , using  results  such  as  those  presented 
in  References  lh , 15,  26,  and  27  as  guides. 

The  concept  of  locking  material  behavior  is  discussed  in  Refer- 
ence 6.  While  the  plastic  locking  strain  or  density  concept  might  be 
a valid  one  for  rocks,  it  is  felt  that  the  actual  locked  plastic  density 
is  not  the  appropriate  value  to  use  for  in  Equation  2.1.  The 

locked  plastic  density  may  be  approached  in  the  material  surrounding 
the  point  of  impact,  but  the  density  increase  will  be  much  less  during 
the  remainder  of  the  penetration  event.  The  pressure-relative  density 
behavior  of  several  rock  types  is  shown  in  Figure  3.6.  For  soils,  the 
"lock-up"  density  or  volumetric  strain  can  be  estimated  as  being  equal 
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to  the  porosity.  However,  due  to  the  strength  of  the  rock  matrix  itself, 
the  porosity  of  rocks  is  not  such  a controlling  factor,  although  it  ob- 
viously affects  the  total  compressibility.  The  porosity  of  rocks  ranges 
from  essentially  zero  to  as  much  as  30  percent.  If  the  impact  pressure, 
average  pressure  in  the  media  during  penetration,  existing  cavity  pres- 
sure, etc.,  can  be  estimated  or  approximated  and  the  initial  density  pq 
is  known,  then  a value  for  pp  can  be  obtained  from  a pressure-relative 
density  relation  (such  as  Figure  3.6). 

The  parameter  e^  has  been  shown  not  to  be  a sensitive  parameter 
in  maximum  penetration  predictions  (Reference  1*).  It  is  postulated  that 
setting  e^  equal  to  zero  for  rocks  is  an  appropriate  assumption. 

3.4  SELECTION  OF  GENERAL  TARGET  MATERIAL  PROPERTIES 

3.4.1  Selection  of  Strength  and  Modulus  Values  for  Representative 
Rock  Groups.  It  was  desired  to  keep  the  parameter  study  in  this  report 
very  general  with  respect  to  the  rock  target  parameters.  Using  the  in- 
tack rock  classification  scheme  presented  in  Section  3.2  as  a guide, 
low-strength  (6000  psi),  medium-strength  (12,000  psi),  and  high-strength 
(24,000  psi)  rock  targets  were  selected,  with  the  strength  values  rep- 
resenting the  median  values  for  Classes  D,  C,  and  B,  respectively 
(Table  3.1).  The  targets  were  characterized  as  being  in  the  average 
modulus  ratio  Class  M with  modulus  ratio  of  350.  Thus  the  three 
targets  would  be  classified  as  BM,  CM,  and  DM.  The  strain-hardening 
modulus  E^  is  set  equal  to  zero  based  on  the  arguments  presented  in 
Section  3.2. 

3.4.2  Selection  of  Initial  Densities  and  Methods  for  Estimating 
the  Compressibilities.  The  remaining  three  parameters  needed  for  the 
three  rock  targets  as  input  to  the  penetration  model  are  the  initial 
densities  pQ  , locked  elastic  volumetric  strains  e^  , and  locked  plas- 
tic volumetric  strains  ep  or  plastic  densities  pp  . While  it  is  not 
true  in  general  that  rocks  with  higher  densities  have  higher  strengths, 
it  is  true  in  many  cases  and  correlations  exist  to  this  effect  (Refer- 
ence 16).  Thus  initial  densities  were  selected  which  were  thought 

to  be  representative  of  the  three  strength  classes  in  an  increasing 
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fashion:  2.1,  2.6,  and  2.9  g/cnr  for  low-,  medium-,  and  high-strength 

targets,  respectively.  The  locked  elastic  volumetric  strain  e^  was 
set  equal  to  zero  for  all  three  rock  targets. 

The  procedure  used  to  estimate  for  the  three  rock  targets  was 

to  estimate  first  the  peak  stress  at  impact  by  graphical  impedance 
2 

matching**  (for  the  maximum  impact  velocity  of  3000  fps  considered  in 
this  report)  in  pressure  versus  particle  velocity  space  (Figure  3.7)  and 
then  the  compressibilities  at  the  peak  stresses  from  the  relations  in 
pressure  versus  PQ/P  space  (Figure  3.8).  The  value  of  for  use  in 

Equation  2.1  is  then  estimated  from  the  compressibility  at  a stress 
level  approximately  one-half  the  peak  stress  for  each  rock  type.  Nei- 
ther the  method  used  here  nor  the  actual  curves  in  Figures  3.7  and  3.8 
are  intended  to  be  rigorous  or  strictly  representative  of  particular 
rocks.  The  curves  are  intended  to  be  representative  of  rocks  in  the 
three  strength  categories  and  are  composites  of  data,  such  as  those 
presented  in  Figure  3.6  and  References  28,  30,  33,  and  3^.  Discussion 
of  the  graphical  impedance  matching  techniques  and  the  equations  con- 
necting the  various  parameters  of  interest  can  be  found  in  many  sources 
(see  References  3^  through  37  for  examples).  Table  3.3  summarizes  the 
assumed  rock  properties  for  DM,  CM,  and  BM  rocks. 


Plane  wave  assumption. 


26 


TABLE  3.1  ENGINEERING  CLASSIFICATION  FOR  INTACT  ROCK 


I.  On 

basis  of  strength 

Uniaxial  Compressive  Strength 

Class 

Description 

psi 

kbar  . (io9 

V cm  / 

A 

Very  high  strength 

Over  32,000 

2.2 

B 

High  strength 

16,000-32,000 

1.1-2. 2 

C 

Medium  strength 

8,000-16,000 

0.55-1.1 

D 

Low  strength 

14,000-8,000 

0.28-0.55 

E 

Very  low  strength 

Less  than  14,000 

0.28 

II.  On 

basis  of  modulus  ratio 

Class 

Description 

a 

Modulus  Ratio 

H 

High  modulus 

ratio 

Over  500 

M 

Average  modulus  ratio 

200-500 

I, 

Low  modulus  ratio 

Less  than  200 

Class!  fy 

■ rock  as  BM,  BH,  BL,  etc. 

n 

Modulus  ratio  = E. cn/o 
tpu  a 

where 

E.c-  = tangent  modulus  at  50 % ultimate 

o - uniaxial  compressive  strength, 
a 

strength ; 

TABLE  3. 

2 ENGINEERING  CLASSIFICATION  FOR  IN  SITU  ROCK 

ROD,  l 

Velocity  Index 

Fracture 
Frequency 
No. /ft 

Description 

0-25 

0-0.2 

6-14.5 

Very  poor 

25-50 

0.2-0. L 

14.5-3 

Poor 

50-75 

0.L-0.6 

3-1.5 

Fair 

75-90 

0.6-0. 8 

1.5-0.75 

Good 

90-100 

0. 8-1.0 

0.75-0 

Excellent 
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TABLE  3.3  ROCK  PROPERTIES 


cd 
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Average  modulus  ratio,  i.e.,  E/Y,  = 350  for  all  three  rock  types. 
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figure  3.1  Engineering  classification  for  intact  rock;  summary  plot  for 
igneous  rocks  (176  specimens,  75  percent  of  points)  (Reference  17). 
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R< ck  strength  based  on  Schmidt  hardness  (Reference  16). 
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Figure  3.1*  Variation  of 
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(Reference  20). 
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Figure  3.6  Presure-relative  density  relations 
for  several  rock  types. 
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CHAPTER  1* 


COMPARISON  OF  PENETRATION  MODEL  PREDICTIONS 
WITH  EXPERIMENTAL  RESULTS 

Ii.l  COMPARISON  01'  PENETRATION  MODEL  PREDICTIONS  USING  GENERAL  ROCK 

PROPERTIES  WITH  EXPERIMENTAL  RESULTS 

In  an  effort  to  gain  insight  concerning  the  reliability  of  the  pen- 
etration model  for  predicting  actual  dynamic  penetration  into  rock,  five 
field  tests  were  modeled  and  penetration  predictions  obtained.  The 
tests  were  selected  from  those  reported  in  Reference  2.  The  author  of 
the  reference  reported  RQD  values  for  the  targets  and  classified  the 
rock  according  to  the  scheme  presented  in  Table  3.1.  Based  on  the 
classification  of  the  targets,  the  general  rock  properties  listed  in 
Table  3.3  were  used  as  presented,  except  that  if  the  target  was  listed 
as  having  a high  modulus  ratio,  the  value  500  was  used  in  conjunction 
with  the  Y value  to  obtain  E ; likewise,  for  low  modulus  ratio  classi- 
fications, a value  of  200  was  used.  Table  L . 1 presents  the  projectile 
characteristics  and  target  properties  used  in  the  penetration  model  and 
compares  the  predicted  and  measured  depths  of  penetration.  The  projec- 
tile type  and  test  numbers  refer  to  identification  numbers  used  in  Ref- 
erence 2.  The  agreement  between  the  predicted  and  measured  penetration 
depths  is  very  good  in  most  cases.  It  is  emphasised  again  that  the  rock 
properties  used  were  those  of  Table  3.3  for  low-,  medium-,  and  high- 
strength  rock  and  that  actual  stress-strain  and  strength  data  were  not 
available.  The  fact  that  adequate  simulation  could  be  obtained  from 
such  limited  rock  property  information  is  fortunate,  since  the  military 
targeteer  rarely  has  detailed  material  property  data  for  specific 
targets  available  to  him. 

The  worst  agreement  was  obtained  for  Test  No.  120-127,  for  which 
the  prediction  is  36  percent  lower  than  the  measured  value  for  penetra- 
tion into  Madera  limestone.  For  this  case  the  RQD  at  the  site  was 
not  reported.  For  Test  No.  120-106,  the  initial  prediction  was  about 
25  percent  too  low.  Based  on  Judgment,  it  was  felt  that  the  low  value 
of  RQD  (32  percent)  for  the  weathered  granite  at  the  test  site  was  due 
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to  vertical  as  well  as  to  horizons  il  fractures  (contrasted  to  a sedimen- 
tary rock,  for  which  a low  value  < f RQD  might  be  due  entirely  t.r  hori- 
zontal fractures).  Since  vertical  fractures  would  obviously  enhance 
penetration  and  reduce  strength,  a second  calculation  was  attempted  with 
the  value  of  Y lowered  to  U000  psi , which  is  the  lower  limit  for  low- 
strength  rocks  in  Table  3.1.  This  attempt  yielded  results  which  agreed 
much  better  with  those  obtained  in  the  experiment  (second  entry  for 
Test  No.  120-106  in  Table  ’4.1).  In  similar  fashion,  a second  calcula- 
tion for  Test  No.  120-127  with  Y lowered  to  8000  psi  (the  lower  limit 
of  CH  rocks)  yielded  much  better  agreement  with  the  experimental  result 
(second  entry  for  Test  No.  120-127). 

For  the  limited  number  of  comparison  cases,  the  good  agreement 
shown  in  Table  U . 1 could  be  due  to  fortuitous  combinations  of  material 
parameters.  Thus  it  is  profitable  to  obtain  upper  and  lower  limit 
penetration  predictions  for  each  case  and  then  compare  these  with  the 
measured  penetration  depths.  This  comparison  is  shown  in  Figure  U.l, 
in  which  predicted  upper  and  lower  limit  depths  of  penetration  are  plot- 
ted versus  measured  depth  of  penetration  (the  perfect  correlation  line 
is  shown  for  reference).  The  upper  and  lower  limits  of  penetration  for 
each  case  were  calculated  using  the  lower  and  upper  limits,  respec- 
tively, of  the  compressive  strength  for  a given  intact  rock  classifica- 
tion (i.e.,  a factor  of  two  variation  in  strength).  Furthermore,  since 
the  compressibility  characteristics  of  the  target  materials  were  not 
known,  it  was  assumed  that  the  upper  limits  of  compressive  strengths 
correspond  to  1 percent  and  the  lower  limits  to  20  percent  compressi- 
bility. Any  other  reasonable  combinations  of  material  properties  will 
result  in  penetration  depths  which  fall  within  these  limits.  As  shown 
in  Figure  U.l,  the  use  of  the  rather  extreme  range  of  properties  still 
produces  ranges  of  predicted  penetration  depths  which  are  acceptable 
for  many  practical  applications. 

f 

U.2  TARGET  PARAMETER  STUDY  OF  SANDIA  TEST  NO.  120-106 

In  order  to  determine  if  the  success  achieved  by  changing  the  yield 
strength  of  the  target  (E  adjusted  to  keep  modulus  ratio  at  500)  and 
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keeping  other  target  properties  constant  was  accidental,  a parameter 
study  of  Sandia  Test  No.  120-106  was  conducted  to  illustrate  the  sensi- 
tivity of  penetration  predictions  to  all  the  target  properties.  The 
projectile  parameters  were  not  varied:  weight  W = 6l3  pounds,  diame- 

ter D = 8 inches,  sectional  pressure  W/A  * 12.2  psi , and  CRH  = 9*25  . 
Figure  k.2  illustrates  the  effect  of  varying  Y within  the  range  1*000- 
8000  psi,  using  E , pq  , and  from  Table  3.3  for  a low-strength 

rock.  The  penetration  data  in  Figure  1*.2  are  replotted  in  Figure  1*.3 
for  selected  velocities  as  a function  of  modulus  ratio,  where  in  this 
case  the  modulus  ratio  varies  due  only  to  Y since  E is  held  con- 
stant. In  Figure  l*.l*,  E is  varied  by  a factor  of  20,  and  the  effect 
on  penetration  is  seen  to  be  very  small  (less  than  5 percent  over  the 
entire  velocity  range ) . Thus , for  a factor  of  two  change  in  Y , the 
effect  on  predicted  penetration  depth  for  a given  velocity  will  be  of 
the  order  of  1*0  percent,  while  E has  a minimal  effect  (at  least  for 
a factor  of  10  change  in  E). 

The  effect  of  rock  target  density  and  compressibility  on  penetra- 
tion is  investigated  in  Figures  1*.5  through  1*.7  for  constant  values  of 
E and  Y . A constant  value  of  the  quantity  Pp/P0  is  maintained  in 
Figure  1*.5  while  varying  both  pQ  and  p^  ; the  effect  on  penetration 
is  seen  to  be  minimal  at  low  velocities  and  approaches  a spread  of  only 
7 percent  about  the  median  penetration  depth  at  1900  fps.  The  effect 
of  varying  the  compressibility  while  keeping  pq  constant  is  shown  in 
Figure  1*.6,  and  Figure  h.'f  shows  the  effect  of  varying  compressibility 
while  keeping  p^  constant.  While  it  is  obviously  important  to  choose 

realistic  values  for  p and  p , it  is  the  ratio  p /p  which  ap- 

o p p o r 

pears  to  have  a more  significant  effect  on  predicted  depth  of  penetra- 
tion (for  a given  Y).  As  shown  in  Figure  1*.6,  in  the  upper  range  of 
compressibility  (for  example,  1.11  £ pp/p0  5 1.35),  a variation  of 
10  percent  in  Pp/P0  affects  the  predicted  depth  of  penetration  by 
the  order  of  10  percent.  For  the  relatively  incompressible  range  (for 
example,  1.00  < Pp/P0  - 1.01),  a variation  of  only  1 percent  affects 
the  predicted  det  \ of  penetration  by  the  order  of  15  percent.  For 
the  range  of  velocity  of  1000-3000  fps,  the  value  of  Pp/P0  should  be 


37 


in  thr  range  1.01  < Pp/PQ  s 1.20  for  moat  rocks.  Considering  Fig- 
ure 3.8,  a reasonable  maximum  variation  of  the  estimate  of  Pp/PQ 
should  be  +5  percent.  For  the  above  range  of  Pp/PQ  » this  variation  in 
Pp/pQ  would  produce  maximum  variations  in  the  predicted  depth  of  pene- 
tration of  the  order  of  +10  percent. 

For  Test  No.  120-106,  the  RQD  in  situ  is  given  as  32  percent; 
using  the  general  parameter  value  (Table  3.2)  for  low-strength  rock  of 
2.1  x 10^  psi  as  the  intact  modulus,  the  in  situ  modulus  can  be  esti- 
mated using  Figure  3.6  (upper  curve)  to  be  0.1*  x 10^  psi.  Reference 
to  Figure  4.4  indicates  that  this  reduction  in  E will  probably  have 
negligible  effect  on  penetration  at  the  impact  velocity  of  860  fps. 

Thus,  it  is  seen  that,  for  this  case,  a reasonable  range  of  variation 
in  Y has  a much  more  significant  influence  than  does  a reasonable  range 
of  variation  of  any  other  parameter  in  the  theory.  It  is  expected  that 
this  same  statement  will  hold  true  for  each  of  the  other  four  rock  pene- 
tration tests  simulated. 

4.3  TENTATIVE  PROCEDURE  FOR  BOUNDING  PENETRATION  DEPTH  IN  ROCK 

It  Is  proposed  that  reasonable  upper  and  lower  bounds  for  penetra- 
tion predictions  into  rock  within  the  range  of  interest  for  EPW's  can 
be  obtained  if  a strength  classification  such  as  "DH"  for  the  intact 
target  material  is  known  or  can  be  inferred  (e.g.,  from  comparison  with 
similar  rocks,  from  aerial  photos,  geologic  maps,  or  literature,  etc.) 
by  using  the  pq  , pp  , and  E values  from  Table  3.3  and  varying  Y 
between  the  limits  indicated  for  the  particular  strength  classification 
in  Figure  3.1.  If  in  addition  the  RQD  is  known  or  can  be  estimated, 
the  E value  can  be  lowered  by  use  of  Figure  3.5  as  previously  shown. 
This  concept  is  illustrated  in  Figure  4.8  for  Test  No.  120-106.  For 
badly  weathered  or  fractured  rock,  this  procedure  may  not  be  valid,  and 
more  work  is  needed  on  correlations  between  RQD  » for  example,  and  in 
situ  strength.  The  measured  penetration  depth  for  Test  No.  120-106  lies 
at  the  upper  bound  of  Figure  4,8  (RQD  = 32  percent).  The  cutoff  for 
the  validity  of  this  method  may  be  in  the  range  of  RQD  = 30-40  percent, 
although  as  discussed  in  Chapter  3t  the  type  of  rock  may  influence  the 


38 


limit  of  validity  of  the  method.  While  the  procedure  proposed  here 
needs  further  verification,  it  is  in  general  agreement  with  all  known 
relevant  data  and  is  based  on  a rational  theory  whose  assumptions  and 
limitations  are  defined. 

U.U  EXTENSION  OF  PENETRATION  MODEL  TO  TREAT  THE  COMPLETE  TARGET 

PRESSURE-DENSITY  CURVE 

In  Chapter  3,  for  simplicity  and  generality,  a single  value  of 
locked  plastic  density  was  estimated  for  each  hypothetical  rock  target 
and  used  in  the  penetration  model  for  the  entire  velocity  range  con- 
sidered. The  effect  of  target  compressibility  on  depth  of  penetration 
predictions  has  been  discussed  previously  in  this  report.  The  effect 
of  compressibility  on  deceleration-time  histories  will  be  discussed  in 
this  section.  A better  way  of  treating  the  compressibility,  which  is 
still  within  the  framework  of  the  constant  locked  plastic  density  con- 
cept, would  be  to  input  a value  for  , estimated  as  in  Section  3.3.2, 

for  example,  for  each  impact  velocity.  The  idea  of  using  a value  for 
compressibility  at  a stress  equal  to  one-half  the  impact  stress  at  the 
maximum  velocity  of  interest  was  arbitrary  and  a representation  of  the 
fact  that  the  stress  in  the  target  media  drops  very  rapidly  following 
impact  and  in  fact  is  probably  much  below  one-half  the  peak  stress 
value  for  most  of  the  penetration  event.  To  improve  the  treatment  of 
compressibility  in  the  penetration  model,  a logical  and  desirable  ex- 
tension is  to  input  the  correct  pp  for  the  existing  cavity  pressure 
at  each  integration  cycle,  i.e.,  the  complete  target  pressure  versus 
density  curve  should  be  an  input  to  the  model.  Since  the  cavity  pres- 
sure is  available  from  the  calculations  in  the  computer  codes,  the 
extension  to  treat  the  complete  target  pressure-density  relation  was  not 
difficult  and  was  carried  out  as  pa/'t  of  this  study. 

U.5  EVALUATION  OF  EFFECT  OF  PRESSURE-DENSITY  CURVE  MODIFICATION 

The  improvements  made  are  expected  to  have  their  greatest  influence 
on  the  deceleration-time  history.  Regrettably,  no  deceleration  records 
were  available  for  the  five  events  simulated  earlier  in  this  chapter. 
However,  Reference  38  presents  calculated  deceleration-time  histories 
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for  particular  points  in  deformable  projectiles  as  part  of  the  results 
of  numerical  solutions  to  projectile  penetration  into  two  rock  types. 
These  numerical  solutions  were  obtained  by  using  a large,  multipurpose, 
two-dimensional,  Lagrangian  finite  difference  computer  code.  The  two 
rock  types  were  welded  tuff  (low-strength  rock)  and  Madera  limestone 
(medium-strength  rock).  The  results  of  these  two  calculations  are  used 
to  test  the  penetration  model  as  presented  in  this  report  and  as  modi- 
fied to  treat  the  complete  target  pressure-density  relation  using  pro- 
jectile and  target  properties  presented  in  Reference  38.  The  target 
hydrostats'' are  modeled  in  Reference  38  using  the  Mie-Grtineisen  equation 
of  state: 

p * „P°  ° ,g  [x  » r ,|(1  - n)j  ♦ VoEi 

(1  - qn)  L -• 

where 

Pq  = initial  density 

c = initial  velocity  of  sound 
o 

n = 1 - PQ/P 

q = material  constant  relating  shock  velocity  to  sound  and 
particle  velocities 

Y = Grtineisen  constant 
o 

E.  = internal  energy 
p = density 

Table  U.2  (data  taken  from  Reference  38)  summarizes  the  target  and 
projectile  properties  for  the  two  penetration  events.  In  the  finite 
difference  code  calculation,  very  large  yield  strengths  were  input  to 
assure  that  no  plastic  deformations  took  place  in  the  projectile.  Body 
forces,  heat  conduction,  and  energy  sources  were  also  neglected.  In 
the  present  penetration  model,  a rigid  projectile  is  assumed  and  the 
other  factors  are  not  considered.  Figures  h.9  and  ^.10  are  pressure 
versus  density  relations  obtained  from  Equation  l».l,  using  the  data  in 
Table  U.2\  Figure  l*.ll  is  a higher  pressure  extension  of  Figure  h.9. 

Figure  h . 12  is  a reproduction  of  a plot  from  Reference  38  showing 
the  finite  difference  code  computed  deceleration  versus  time  response 
for  the  welded  tuff  penetration  calculation  of  a particle  located  at 
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the  tip  of  the  projectile  nose.  Superimposed  on  this  is  the  rigid 
body  prediction  from  the  penetration  model  calculation  (dashed  line). 

In  Figure  1».13  the  penetration  model  deceleration-,  velocity-,  and 
displacement-time  histories  (dashed  lines)  are  superimposed  on  the  time 
histories  calculated  by  the  finite  difference  code  for  a particle  lo- 
cated on  the  centerline  of  the  projectile  at  a distance  equal  to  the 
nose  length  frou  the  tip.  As  can  be  seen  in  these  figures,  the  penetra- 
tion model  deceleration -time  histories  are  quite  similar  to  the  lower 
frequency  portion  obtained  from  the  deformable  body  finite  difference 
calculations.  These  high-frequency  oscillations  represent  particle 
accelerations  and  decelerations  due  to  stress  wave  propagation  in  the 
projectile,  which  of  course  cannot  be  simulated  with  a rigid  body  pene- 
tration model.  As  shown  in  Figure  1<.13,  the  effect  of  differences  in 
the  acceleration  curves  has  a negligible  effect  on  velocity-  and 
displacement-time  histories  from  the  two  calculational  methods.  It  is 
interesting  to  note  that  the  observed  final  penetration  depth  for  the 
test  simulated  in  Figure  U.13  was  k.95  feet.  This  value  is  in  close 
agreement  with  the  calculated  results. 

As  discussed  in  Reference  38,  although  the  peak  decelerations  in 
Figures  U.12  and  U.13  for  the  finite  difference  calculation  are  drasti- 
cally different,  the  decelerations  for  both  curves  should  oscillate 
about  the  common  "rigid  body  deceleration"  of  the  projectile;  and  it  is 
evident  that  the  agreement  of  the  penetration  model  prediction  of  the 
rigid  body  motion  with  that  predicted  from  the  finite  difference  code 
results  is  excellent  in  terms  of  general  shapes  of  the  time  histories, 
magnitudes,  and  duration  of  the  event. 

Figure  U.lU  presents  the  time  histories  for  the  limestone  pene- 
tration calculated  with  the  penetration  model  superimposed  on  the 
finite  difference  code  generated  time  histories  (base  figure  reproduced 
from  Reference  38).  For  this  case,  the  time  histories  are  for  a par- 
ticle on  the  projectile  centerline  at  a distance  of  2 feet  from  the 
nose  tip.  Again  the  penetration  model  predictions  are  in  excellent 
agreement  with  the  finite  difference  results.  The  maximum  depth  of 
penetration  observed  in  the  event  simulated  here  was  3.6  feet.  Again 
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both  calculations  are  in  reasonably  good  agreement  with  this  value. 

Also,  both  the  penetration  +heory  and  th  ' finite  difference  deceleration- 
time results  are  in  qualitative  agreement  with  the  experimental 
deceleration-time  history  shown  in  Figure  U . ll»  (this  represents  the 
only  available  complete  deceleration-time  record  for  projectile  penetra- 
tion into  a nonlayered  rock  target). 

To  illustrate  the  importance  of  the  details  of  target  compressibil- 
ity on  the  prediction  of  deceleration-time  histories,  a parameter  study 
was  conducted.  The  results  are  shown  in  Figures  ^.15  and  U.16.  Fig- 
ure k.15  shows  the  effect  on  penetration  versus  time  plots  of  using  the 
constant  locked  plastic  density  model  with  various  compressibility 
estimates  for  the  welded  tuff  case.  The  plot  for  0.5  percent  compress- 
ibility agrees  approximately  in  final  penetration  depth  and  event  dura- 
tion with  the  prediction  based  on  the  complete  pressure-density  curve 
as  shown  in  Figure  14.13.  Figure  I4.16  shows  deceleration  versus  time 
plots  for  the  5.3  and  the  0.5  percent  compressibilities  along  with  the 
deceleration-time  curve  from  Figure  * .13.  The  concept  begins  to  emerge 
that  if  the  "right"  constant  locked  /plastic  density  value  could  be 
selected,  then  the  simple  constant  locked  plastic  density  model  could  be 
used  to  obtain  results  in  good  agreement  with  those  obtained  when  the 
complete  pressure-density  relation^is  known  and  used. 

If  the  details  of  the  deceleration-time  history  of  rock  penetration 
events  are  of  interest,  it  is  clearly  important  to  choose  the  "right" 
pp  . How  should  one  go  about  estimating  a "right"  constant  locked  plas- 
tic density  for  a particular  penetration  problem?  Another  estimate  of 
target  compressibility  can  be  obtained  by  considering  the  contribution 
to  the  dynamic  cavity  pressure  expressed  by  the  term  fCeJp^B^V  , 
where  V is  the  velocity  (see  Equation  2.1).  By  using  the  impact  veloc- 
ity Vq  in  this  term,  an  estimate  is  obtained  of  the  peak  dynamic  cav- 
ity pressure  contribution  from  this  term.  For  impact  and  penetration 
into  the  medium-strength  rock  of  Chapter  3 at  a velocity  of  1000  fps, 
with  CRH  = 6 , W = 750  pounds,  and  D = 8 inches,  the  calculated 
peak  dynamic  cavity  pressure  is  about  2 kbar.  Following  the  procedure 
of  estimating  the  compressibility  at  half  of  this  peak  stress,  which  was 
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proposed  in  connection  with  Figure  3.8,  yields  approximately  1 percent 
compressibility  compared  to  (l)  the  5 percent  compressibility  estimated 
from  the  same  procedure  but  using  peak  stresses  at  impact  estimated 
from  impedance  matching  for  Vq  * 1000  fps  and  (2)  the  13  percent  com- 
pressibility used  in  Table  3.2  (based  on  impedance  matching  for  the 
upper  limit  of  VQ  = 3000  fps).  Figure  1+.17  shows  the  time  histories 
for  this  case  for  a 1 percent  compressibility  and  for  the  case  in  which 
the  entire  hydrostat  of  Figure  3.8  was  used.  The  event  duration,  depth 
of  penetration,  and  peak  deceleration  agree  very  closely.  The  general 
effect  of  using  the  model  with  constant  locked  plastic  density  with  a 
compressibility  which  is  too  large  is  to  overestimate  the  depth  of  pene- 
tration and  the  event  duration  and  to  underestimate  the  peak  decelera- 
tion. A compressibility  of  0.5  percent  (shown  in  Figure  4.16)  is  too 
small,  at  least  for  the  early  part  of  the  event,  and  the  peak  decelera- 
tion is  slightly  overpredicted  compared  with  the  plot  for  the  complete 
pressure-density  relation  calculation. 

It  is  emphasized  that,  if  only  the  final  depth  of  penetration  is  of 
interest,  then  the  need  to  select  the  "right"  pp  is  net  of  overriding 
importance,  and  a simple  method  such  as  the  impedance  matching  procedure 
will  be  adequate.  Based  on  these  considerations,  the  original  rock 
property  values  given  in  Table  3.3  were  used  in  the  parameter  study  in 
Chapter  5,  since  penetration  depth  was  the  primary  dependent  variable 
of  interest.  However,  the  penetration  predictions  to  be  presented  in 
Chapter  5 should  be  viewed  as  upper  bound  estimates,  particularly  at 
the  lower  impact  velocities,  since  the  compressibilities  are  somewhat 


TABLE  It.l  COMPARISON  OF  PENETRATION  MODEL  PREDICTIONS  WITH  FIE1JD  PENETRATION  TEST  RESULTS 
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TABLE  U.2  PROPERTIES  FOR  PROJECTILE  PENETRATION 
INTO  TUFF  AND  LIMESTONE 


Welded  Tuff  Penet rat Ion 


Tuff  Properties 


pQ  * 1.85  g/cm  * 115  pcf 


Poisson's  ratio  * 0.1 

= 3.8 
= 5510  psi 


Y = 3.8  x 10^  dynes/cm2 


5 


* 2.05  * 10  cm/ sec 
= 6725  fps 


Y * 2.00 
o 


1.1*9 


Pro.jectile  Properties 
D * 22.86  cm  * 9 in 
W = 1*.55  * 105  g = 1000  lb 
W/A  = 15.76  psi 


Nose 

CRH 


V = 21,181*  cm/sec  * 695  fps 
o 


Madera  Limestone  Penetration 


Limestone  Properties 


p =2.69  g/czn  = 168  pcf 
o 1 


Poisson's  ratio  = 0.32 


Y = 9.1*1*  x 10®  dynes/cm2 
= 13,690  psi 

5 

: = 3.1*  x 10  cm/sec 

0 « 11,150  fps 

r * 2.00 

q s 1.27 


Projectile  Properties 

D = 20.32  cm  = 8 in 

W = 3.055  x 105  g 
* 673.6  lb 

W/A  -13.1*  psi 


Nose 

CRH 


9.25 


Vq  = 17,371*  cm/sec  * 570  fps 


PREDICTED  DEPTH  OF  PENETRATION,  IN 


I 


Figure  k.l  Measured  versus  predicted  depth  of 
penetration  into  rock  targets. 


1+6 


PENETRATION,  IN  PENETRATION 


500  TOO  *00  1100  1300  1500  1700  1500  2100 

IMPACT  VELOCITY,  TP S 


Figure  4.2  Effect  of  varying  Y on  penetration  versus 
impact  velocity  plot  (Test  No.  120-106). 
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Figure  4.3  Penetration  versus  modulus  ratio  for  constant  E. 
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Figure  U . L Effect  of  modulus  on  penetration. 
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Figure  1+.5  Effect  of  pQ  and  p on  penetration  while 

maintaining  p /p  constant.  ^ 
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Figure  h.6  Effect  of  p /p  on  penetration  while  maintaining 
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o 


MOO  1)00  1500 

IMPACT  VELOCITY,  TPS 


Figure  U.7  Effect  of  p /p  on  penetration  while  keeping 
. , p o 

p constant. 
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Figure  U.10  Pressure-density  relation  for  Madera 
limestone  based  on  Equation  b.l. 
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Figure  14.12  WES  rigid  body  acceleration  versus 
time  prediction  superimposed  on  the  Sandia 
finite  difference  results  for  particle  located 
on  nose  tip  of  projectile;  tuff  penetration. 
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Figure  U.13  WES  rigid  body  model  time  histories 
superimposed  on  the  finite  difference  time 
histories  for  particle  on  centerline  located  one 
nose  length  from  nose  tip;  tuff  penetration. 
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Figure  U.lt  WES  rigid  body  model  time  histories  superimposed 
on  the  finite  difference  time  histories  for  a particle  on 
centerline  located  2 feet  from  nose  tip;  limestone  penetration 
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Figure  4.15  Constant  locked  plastic  density  model  predictions 
for  various  compressibilities;  tuff  penetration. 


TIME,  MSEC 

Figure  U.l6  Comparison  of  complete  target  pressure-density 
model  with  two  constant  locked  plastic  density  model  pre- 
dictions of  deceleration  versus  time. 
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CHAPTER  5 


GENERAL  ROCK  PENETRATION  PARAMETER  STUDY 

5.1  SYSTEMATICS  OF  PARAMETER  STUDY 

The  rock  properties  used  in  the  rock  penetration  parameter  study 
to  be  discussed  in  this  chapter  were  summarized  in  Table  3.3.  In  the 
remainder  of  this  report,  the  three  rock  parameter  groups  will  be  in- 
dicated by  L (low  strength),  M (medium  strength),  and  H (high 
strength).  Table  5.1  summarizes  the  ranges  of  projectile  parameters 
to  be  considered.  Obviously,  it  was  unnecessary  and  impractical  to 
consider  all  possible  combinations  of  rock  and  projectile  parameters; 
only  those  combinations  of  projectile  weight  and  diameter  which  give 
values  of  W/A  in  the  range  of  10-15  psi  are  considered,  as  this  is 
probably  the  reasonable  range  to  be  expected  for  EPW's.  Table  5.2 
summarizes  some  of  the  combinations  considered. 

The  systematics  followed  during  the  general  parameter  study  are  out- 
lined in  Table  5.3.  In  each  case,  the  final  depth  of  penetration  was 
calculated  using  the  penetration  model  presented  in  Chapter  2 for  ini- 
tial velocities  ranging  from  1000  to  3000  fps  in  increments  of  250  fps. 

It  must  be  emphasized  that  the  projectiles  are  considered  as  rigid 
bodies,  and  hence  the  range  of  velocities  may  not  seem  realistic,  since 
it  is  known  that  structural  failure  of  the  projectile  may  occur  at  ve- 
locities considerably  below  3000  fps  and  that  the  critical  velocity  for 
failure  depends,  of  course,  on  target  strength.  However,  as  the  state 
of  the  art  in  projectile  design  advances  and  as  more  exotic  materials 
become  available  for  use,  the  upper  bound  of  survivable  velocities  for 
impact  and  penetration  into  rock  may  reasonably  be  expected  to  increase. 
The  next  three  sections  present  the  results  of  the  parameter  study,  with 
a discussion  of  the  results  given  in  Section  5.5. 

5.2  TIME  HISTORIES 

Since  the  penetration  equation  presented  in  Chapter  2 (Equation  2.1) 
is  a differential  equation,  the  model  can  be  used  to  produce  time 


57 


histories  of  penetration,  velocity,  and  deceleration1  for  the  projectiles 
as  well  as  final  depth  of  penetration.  In  order  to  illustrate  the  kine- 
matics predicted  by  the  model,  typical  complete  time  histories  are  pre- 
sented for  Cases  IB1  and  IIIB3  at  an  impact  velocity  of  2000  fps  in  Fig- 
ures 5.1a  and  5.2,  respectively.  Fig  are  5.1b  compares  the  deceleration 
versus  time  plots  for  the  three  rock  strengths  of  Cases  IB1  through  IB3. 

5.3  PENETRATION  VERSUS  IMPACT  VELOCIT'  PLOTS 

Figures  5.3  through  5.1^  are  final  depth  of  penetration  versus  im- 
pact velocity  plots  for  Cases  I through  IV  of  Table  5.3;  each  figure 
contains  plots  for  the  three  rock  strength  groups  with  other  parameters 
held  constant. 

5.*4  ADDITIONAL  PARAMETER  STUDY  PLOTS 

Using  the  data  from  Figures  5.3  through  5.1^  and  the  results  of 
other  cases  as  required,  further  plots  were  made  illustrating  the  ef- 
fects of  various  parameters  on  predicted  penetration.  Figure  5-15  il- 
lustrates the  effect  of  projectile  diameter  and  hence  sectional  pressure 
W/A  on  the  penetration  versus  impact  velocity  plot  for  medium-strength 
rock  and  a given  CRH  and  W , i.e.,  W/A  varies  due  only  to  a varying 
projectile  diameter.  The  effect  of  CRH  on  penetration  is  shown  in 
Figures  5.16a  and  5.16b  for  low-  and  medium- strength  rock,  respectively, 
and  for  a given  projectile  W/A  at  three  impact  velocities.  Fig- 
ure 5.16c  shows  the  effect  of  CRH  on  predicted  penetration  for  the 
three  rock  types  at  a given  impact  velocity.  Figures  5.17a  through  5.17f 
all  show  the  effect  of  W/A  on  predicted  penetration  for  various  com- 
binations of  impact  velocity  and  rock  type. 

5.5  DISCUSSION  OF  RESULTS 

5.5.1  Time  Histories.  The  time  histories  in  Figures  5.1  and  5.2 


The  term  deceleration  is  used  to  refer  to  a negative  acceleration. 
Since  the  projectiles  are  treated  as  rigid  bodies  (i.e.,  no  wave 
motions  considered  in  the  projectile)  with  the  positive  z-axis  down- 
ward, all  accelerations  will  be  negative. 
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illustrate  the  kinematics  of  the  penetration  model  using  a constant 
value  for  the  locked  plastic  density.  Figure  5.1b  illustrates  the  con- 
cept of  increasing  peak  deceleration  and  decreasing  duration  of  the 
penetration  event  with  increasing  rock  strength  such  that  the  area  under 
the  deceleration  versus  tim_  plots  remains  constant  for  a given  impact 
velocity.  It  is  important  to  realize  that  although  the  penetration 
predictions  given  by  using  the  constant  value  for  locked  plastic  density 
have  been  shown  to  agree  closely  with  experimental  results  (see  Refer- 
ences U,  6,  7.  and  9 and  also  Chapter  1*  of  this  report),  the  decelera- 
tion versus  time  predictions  (rigid  body)  do  not  in  many  cases  reproduce 
experimentally  measured  records.  Deceleration  versus  time  records  from 
actual  penetration  tests  generally  do  not  exhibit  the  sharp  acceleration 
peak  followed  by  a rapid  drop  to  a much  lower  value  which  is  typical 
of  the  deceleration-time  curves  calculated  by  the  penetration  model 
using  a constant  . Rather  the  experimental  records  show  a rise  to 
a peak  (as  the  nose  imbeds)  followed  by  a much  smaller  decrease  during 
the  remainder  of  the  penetration  event.  Indeed,  the  experimental  rec- 
ords in  many  cases  resemble  a step  pulse  of  magnitude  equal  to  the 
average  deceleration  A during  the  event  (determined  from  energy-work 
equality)  of  duration  t ; this  is  shown  as  a dashed  line  in  Figure  5.2, 
where 


A 


t = 


(5.1) 


and  Vq  , g , and  P are  the  impact  velocity,  acceleration  of  gravity, 
and  maximum  depth  of  penetration,  respectively.  In  cases  where  the 
actual  deceleration  versus  time  history  can  be  approximated  by  a step 
pulse,  the  peak  and  average  decelerations  may  not  differ  by  more  them 
a factor  of  two  (Reference  7).  Deceleration  versus  time  predictions 
were  discussed  in  Chapter  U,  where  it  was  shown  that  using  a realistic 
nonlinear  hydrostat  instead  of  a constant  locked  plastic  density  pro- 
vides more  realistic  deceleration-time  histories.  However,  since  this 
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parameter  study  was  conceived  primarily  to  study  the  effects  of  inde- 
pendent variables  on  the  maximum  depth  of  penetration,  and  since  P is 
not  seriously  changed  by  using  a constant  pp  instead  of  the  actual 
hydrostat,  the  simple  approach  used  in  this  study  is  adequate. 

5.5.2  Penetration  Versus  Impact  Velocity  Plots.  The  first  and 
perhaps  most  obvious  comment  regarding  Figures  5. 3-5.1**  is  that  for  a 
given  velocity,  penetration  decreases  with  increasing  target  strength 
and  modulus  (maintaining  E/Y  = 350  constant).  For  all  three  target 
strengths,  penetration  increases  with  increasing  velocity  as  expected, 
although  the  curves  tend  to  flatten  with  increasing  velocity.  The 
curves  for  the  three  rock  targets  are  seen  not  to  be  parallel  for  any 
of  the  cases  considered,  i.e.,  the  numerical  difference  between  the  pen- 
etration into  low-strength  rock  P^  and  into  medium-strength  rock  P^ 
does  not  remain  constant  with  increasing  velocity,  and  the  same  is  true 
for  the  difference  between  the  medium- strength  P^  and  high-strength 

P„  penetration  versus  velocity  curves.  Nor  does  the  percentage  de- 
H 

crease  in  penetration  remain  constant,  but  is  seen  to  be  a slowly  de- 
creasing function  of  velocity. 

Figure  5.18  is  a summary  plot  showing  the  percent  decrease  in  pene- 
tration versus  impact  velocity  for  all  the  caseB  considered  in  this 
report  (CRH  range  2-10  and  W/A  range  10-15  psi).  This  figure  is  sig- 
nificant in  that  it  illustrates  qualitatively  the  effect  of  target 
strength  on  final  depth  of  penetration  predictions.  Thus  the  relative, 
qualitative  effect  shown  in  Figure  5.18  is  reasonable  and  may  be  more 
believable  than  the  quantitative  depth  predictions.  Three  well-defined 
bands  are  identified,  with  the  upper  band  representing  the  percent  de- 
crease in  penetration  between  the  low-  and  high-strength  rocks 
[(PL  * V/PJ  x 100  , the  middle  band  representing  the  percent  decrease 
in  penetration  between  the  low-  and  medium-strength  rocks 
[(pt,  - V/pl]  x 100  , and  the  lower  band  representing  the  percent  de- 
crease in  penetration  between  the  medium-  and  high-strength  rocks 
£(Ph  - * 100  • Figure  5.18  illustrates  in  summary  fashion  that 

the  relative  effect  of  increasing  rock  strength  on  decreasing  the 
depth  of  penetration  at  a given  velocity  is  not  strongly  sensitive  to 
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variations  in  CRH  or  W/A  . However,  this  is  only  a relative  effect 
and  does  not  imply  that  for  a given  rock  strength,  penetration  depth  is 
unaffected  by  CRH  or  W/A  . As  an  example  of  the  practical  value  of 
this  information,  it  is  only  necessary  to  have  an  experimental  penetra- 
tion depth  for  one  rock  target  for  a given  projectile  in  order  to  esti- 
mate penetration  depths  for  that  projectile  into  other  rock  targets. 

5.5.3  Effects  of  W/A  and  CRH  . Figures  5.15-5.17  examine  in 
detail  the  effects  of  CRH  and  W/A  on  final  depth  of  penetration  for 
varying  impact  velocities  and  rock  strengths.  The  effect  of  varying  the 
projectile  diameter  for  penetration  into  a medium-strength  rock  while 
maintaining  W and  CRH  constant  is  illustrated  in  Figure  5.15;  W/A 
varies  from  10  to  17.7  psi.  The  increase  in  penetration  over  the  veloc- 
ity range  considered  is  about  70  percent  for  a decrease  in  diameter  from 
8 to  6 inches.  Thus,  it  is  seen  that  for  a W/A  of  10  psi  at  an  impact 
velocity  of  1000  fps,  increasing  the  W/A  from  10  to  17.7  psi  produces 
approximately  the  same  increase  in  penetration  as  does  increasing  the 
impact  velocity  from  1000  to  1500  fps  and  maintaining  W/A  at  10  psi. 
Considerations  such  as  this  should  be  of  great  value  in  providing  guid- 
ance for  penetrator  designers. 

The  effect  of  varying  CRH  is  illustrated  in  Figure  5.16  for  the 
three  rock  targets  and  at  three  impact  velocities.  CRH  is  seen  to 
have  little  effect  on  depth  of  penetration  for  a given  impact  velocity 
and  rock  strength;  the  effect  is  greatest  for  high  impact  velocity  and 
low  rock  strength.  The  curve  for  Vq  * 1000  fps  in  Figures  5.16a  and 
5.16b  is  seen  to  be  virtually  flat  over  the  range  of  CRH  values  (2-10). 
Likewise,  the  curve  in  Figure  5.16c  for  high-strength  rock  at 
Vq  = 2000  fps  shows  virtually  no  change  in  penetration  with  variation 
in  CRH  . Thus,  within  the  limits  of  parameters  considered  in  this 
study,  the  penetration  model  indicates  that  the  choice  of  projectile 
nose  shape  can  be  based  solely  on  structural  integrity  and  terradynamic 
stability  considerations  and  not  desired  depth  of  penetration  considera- 
tions. The  slight  coupling  between  CRH  and  impact  velocity  and  rock 
strength  does  not  appear  to  be  significant. 
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Figure  5.17  demonstrates  the  effect  of  W/A  on  depth  of  penetra- 
tion of  a projectile  with  a given  nose  shape  (CRH  * 6)  into  the  three 
rock  targets  at  three  impact  velocities.  For  a given  rock  strength, 
sensitivity  of  depth  of  penetration  to  W/A  is  seen  to  increase  as  im- 
pact velocity  increases,  and  for  a given  impact  velocity,  the  sensitiv- 
ity to  W/A  decreases  as  rock  strength  increases.  It  is  interesting 
to  note  that  the  penetration  versus  W/A  relations  shown  in  Fig- 
ures 5.17a  through  5.17c  for  all  three  rock  strengths  appear  to  be 
linear  and  pass  between  the  two  predicted  penetration  depths  for 
W/A  * 12.7  psi  (the  two  points  for  each  rock  strength  are  for  projec- 
tiles with  the  same  W/A  but  different  weights  and  diameters).  However, 
these  points  are  separated  from  each  other  by  less  than  10  inches 
(3-20  percent)  for  all  rock  strengths  and  impact  velocities;  thus  the 
effect  is  small  enough  to  ignore  in  most  cases. 

The  effect  of  projectile  W/A  on  penetration  varies  somewhat  with 
impact  velocity  and  rock  target  being  considered;  however,  for  the  range 
of  impact  velocities  and  rock  strengths  considered,  the  penetration 
model  is  relatively  insensitive  to  variations  in  projectile  CRH  . For 
a medium-strength  rock  at  impact  velocities  of  1000,  2000,  and  3000  fps 
and  CRH  * 6 , the  percent  increase  in  penetration  produced  by  a change 
in  W/A  from  10  to  15  psi  is  1*3,  1*6,  and  50  percent,  respectively. 
Changing  CRH  from  2 to  10  for  penetration  into  a medium-strength  rock 
produces  an  average  increase  of  only  about  16  percent  for  the  three 
impact  velocities. 
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TABLE  5.2  W/A  COMBINATIONS 


W/A 

No. 

W,  lb 

D.  in 

Area 

in2 

W/A 

psi 

1 

250 

5 

ll*.6U 

12.73 

2 

500 

7.5 

UU.18 

11.32 

3 

500 

7 

38. 1*8 

12.99 

1* 

500 

8 

50.27 

9.95 

5 

750 

9 

6U. 51 

11.63 

6 

750 

8 

50.27 

ll*.92 

7 

1000 

10 

78.51* 

12.73 
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TABLE  5.3  PARAMETER  STUDY  SYSTEMATICS 


Case 

Case 

Identification 

Rock 

Identification 

Rock 

No. 

W/A  No. 

CRH 

Type 

No. 

W/A  No.  CRH 

Type 

IAI 

6 

2 

L 

IIIB1 

1 6 

L 

2 

M 

2 

M 

3 

H 

3 

H 

IB1 

6 

6 

L 

IIIC1 

1 10 

L 

2 

M 

2 

M 

3 

H 

3 

H 

IC1 

6 

10 

L 

IVA1 

2 2 

L 

2 

M 

2 

M 

3 

H 

3 

H 

I IAI 

7 

2 

l ; 

IVB1 

2 6 

L 

2 

M 

2 

M 

3 

H 

3 

H 

IIB1 

7 

6 

L 

IVC1 

2 10 

L 

2 

M 

2 

M 

3 

H 

3 

H 

IIC1 

7 

10 

L 

VA1 

: 6 

M 

2 

M 

2 

h 

M 

3 

H 

3 

5 

M 

U 

D « 6.5  in. 

L 

IIIA1 

1 

2 

L 

W * 500  lb 

2 

M 

5 

D * 6 in, 

L 

3 

H 

W = 500  lb 

VIA1 

2 U 

M 

2 

8 

M 

3 

k 

L 

1* 

8 

L 

6h 
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Figure  5.2 


Time  histories  for  Case  IIIB3. 
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Figure  5.3  Case  IA. 


I 


66 


PCNEr«»TlON,  IN  PC  NC  TPATlON 


IMPACT  VCLOCirv,  F PS 

Figure  5.6  Case  IIA. 
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Figure  5.7  Case  IIB. 
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Figure  5.8  Case  IIC. 


Figure  5*9  Case  IIIA 
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Figure  5.15  Effect  of  W/A  or  penetration  versus  impact 
velocity  plots  for  given  rock  type,  CRH  , and  weight. 
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a.  Low-strength  rock:  Vq  * 1000,  2000,  and  3000  fps. 


b.  Medium-strength  rock;  V = 1000,  2000, 
and  3000  fps. 

Figure  5.16  Effect  of  CRH  on  predicted  penetration 
(sheet  1 of  2). 
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Figure  5.18  Percent  decrease  in  penetration  versus  impact 
velocity  for  all  cases  considered  in  this  report  between 
the  low-  and  high-strength  rock,  between  the  low-  and 
medium- strength  rock,  and  between  the  medium-  and  high- 
strength  rock. 


CHAPTER  6 


SUMMARY  AND  CONCLUSIONS 

Projectile  penetration  into  earth  media  is  an  extremely  complicated 
phenomenon,  and  the  state  of  the  art  in  penetration  theory  is  such  that 
the  problem  cannot  be  solved  in  all  its  complexity.  Thompson  and  Bryant 
(Reference  39)  discuss  the  philosophical  and  physical  formulation  of  the 
general  penetration  problem,  and  Chapter  3 of  this  report  presents  an 
overview  of  the  problems  involved  in  characterizing  rock  as  the  target 
medium  for  penetration.  Based  on  a theoretical  penetration  model  dis- 
cussed in  Chapter  2,  which  requires  only  target  material  properties  that 
can  be  determined  from  standard  laboratory  tests,  a study  of  projectile 
penetration  into  rock  is  presented  in  Chapters  ^ and  5. 

Three  rock  targets  are  considered,  with  material  properties  chosen 
to  be  representative  of  general  classes  of  rocks  characterized  as  being 
low-,  medium-,  and  high-strength.  The  three  rock  targets  are  charac- 
terized by  only  four  parameters:  (l)  initial  density,  (2)  locked  plas- 

tic density,  which  is  based  on  estimate  of  compressibility,  (3)  modulus, 
and  (M  strength.  For  target  sites  about  which  only  a minimum  amount 
of  information  is  available,  it  is  proposed  that  reliable  penetration 
depth  estimates  can  be  obtained  if  it  is  possible  (on  the  basis  of 
laboratory  examination  of  cores,  drillers'  logs,  air-photo  interpreta- 
tion, or  geologic  maps  and  literature)  to  classify  the  intact  and  in 
situ  material  according  to  established  engineering  rock  classifications 
(Tables  3.1  and  3.2  and  Figures  3.1  through  3.3,  for  example). 

The  modified  Ross-Hanagud  penetration  theory  (Chapter  2)  with  the 
WES  "recipe"  for  selecting  material  properties  for  rock  (Chapter  3) 
adequately  predicts  depth  of  projectile  penetration  into  rock  (Chap- 
ter U),  When  the  model  is  extended  to  include  realistic  pressure- 
density  relations  for  the  targets,  it  also  appears  to  adequately  simu- 
late rigid  body  deceleration-time  histories  for  relatively  low  impact  ve- 
locities and  low-  to  medium-strength  rocks  (Chapter  U).  Therefore,  the 
conclusions  derived  from  parametric  studies  using  this  model  (Chapters  14 
and  5)  have  relevance  for  the  real  world  problem  of  rock  penetration. 


79 


Using  the  general  rock  properties  proposed  in  Chapter  3,  the  results 
of  five  actual  field  penetration  tests  into  rock  (targets  classified  ac- 
cording to  Tables  3.1  and  3.2)  were  compared  with  the  penetration  model 
predictions  in  Table  U,l,  and  agreement  is  shown  to  be  good.  A parameter 
study  of  one  of  the  field  tests  resulted  in  the  following  conclusions: 

1.  For  a given  value  of  modulus,  variation  of  strength  within  the 
bounds  of  the  classification  categories  (low-strength  "D"  in  this  case) 
can  produce  variations  in  penetration  for  a given  velocity  of  the  order 
of  30-Uo  percent  decrease  in  penetration  for  a factor  of  two  increase 
in  strength. 

2.  For  a given  strength,  variation  in  modulus  by  a factor  of  20 
produces  a negligible  change  in  penetration  over  the  velocity  range  500- 
1900  fps. 

3.  For  realistic  values  of  pQ  and  , variation  of  10  percent 
in  the  estimate  of  Pp/P0  for  the  target  will  produce  a change  in  the 
penetration  depth  prediction  of  the  order  of  10  percent. 

Also,  in  Chapter  U the  penetration  model  is  extended  to  treat  the 
complete  target  pressure-density  relation,  and  good  agreement  is  ob- 
tained with  motion-time  histories  from  two  finite  difference  code  calcu- 
lations of  rock  penetration.  For  these  two  cases,  involving  low-  and 
medium-strength  rock  and  rather  low  impact  velocities  (570  and  695  fps), 
inclusion  of  the  detailed  target  pressure-density  relation  in  the 
theoretical  penetration  model  produces  not  only  good  penetration  depth 
predictions  but  also  decel>-  ation  versus  time  histories  in  general  quali- 
tative and  quantitative  agreement  with  experimental  results.  If  this 
trend  continues  when  compared  with  future  test  data  from  higher  impact 
velocity  field  events,  this  will  prove  to  be  a significant  advancement 
in  the  state  of  the  art  of  rock  penetration  predictions.  Also,  if  time 
histories  from  the  penetration  model  are  improved  by  use  of  the  complete 
pressure-density  curves,  better  results  can  be  anticipated  for  predict- 
ing penetration  into  layered  targets. 

On  the  basis  of  the  parameter  study  of  Chapter  5,  the  following 
comments  and  conclusions  can  be  made: 

1.  Peak  deceleration  increases  with  increasing  rock  strength. 
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2.  Event  duration  decreases  with  increasing  rock  strength. 

3.  Final  penetration  depth  decreases  with  increasing  rock  strength 
at  a given  velocity. 

h.  For  a given  rock  target  strength,  penetration  increases  with 
increasing  impact  velocity,  although  the  curves  tend  to  flatten  with 
increasing  velocity. 

5.  For  a given  velocity  and  rock  strength,  penetration  increases 
only  slightly  with  projectile  CRH  , and  the  effect  is  greatest  for 
higher  velocities  and  low  rock  strengths.  Thus,  the  theory  suggests 
that  for  rock  penetration,  projectile  nose  shape  (for  ogives)  can  be 
based  solely  on  consideration  of  structural  integrity  and  stability, 
and  not  desired  depth  of  penetration  considerations. 

6.  The  effect  of  sectional  pressure  W/A  on  penetration,  over 
the  range  considered,  is  more  significant  than  CRH  . 

a.  Penetration  increases  with  increasing  W/A  . 

b.  As  with  CRH  , the  effect  is  more  significant  for  higher 
velocities  and  low  rock  strengths. 

c.  The  information  contained  in  Figures  5«15  and  5.17  should 
be  of  value  to  designers  in  determining  practicality  of  achieving  in- 
creased penetration  depths  by  increased  impact  velocity  and/or  W/A  . 

Further  work  is  needed  in  the  following  areas  and  is  planned  for 
future  study  under  this  project. 

1.  Effect  of  in  situ  rock  condition  (RQD  , for  example)  on  mate- 
rial properties  for  input  to  the  penetration  model. 

2.  Comparison  of  penetration  model  time  histories  with  finite 
difference  calculations  and  new  experimental  results  of  penetration 
into  rock  as  they  become  available. 

3.  Comparison  of  constant  plastic  density  calculations  with  com- 
plete hydros  tat  calculations  of  penetration  into  rock  in  order  to  de- 
velop better  methods  of  estimating  compressibilities  in  cases  where  the 
hydrostat  is  not  known. 

1.  Oblique  impact  and  consideration  of  angle  of  attack. 

5.  Inclusion  of  projectile  material  properties  in  the  penetration 
model,  if  possible  without  greatly  complicating  the  model. 
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APPENDIX  A 


NOTATION 


A 

A 

BM 


c 

o 


CM 

CRH 


Deceleration 
Average  deceleration 

Classification  of  low-strength  rock  targets 
1 - 61/3 
3/2  - (1  + 

Initial  velocity  of  sound 

Classification  of  medium- strength  rock  targets 
Caliber  Radius  Head 


D 

DM 

e . 

l 


e 

P 


E 


E. 

l 


E. 

msitu 


intact 

E, 


E 


v;o 

EPW 


Diameter  of  projectile 

Classification  of  high-strength  rock  targets 

Volumetric  strain  related  to  the  jump  from  zero  to  finite 
stress  (elastic  region)  in  the  stress-strain  curves  in 
Figure  2.1 

Volumetric  strain  related  to  the  jump  from  the  elastic  to 
the  plastic  region  of  the  uniaxial  stress  (i.e.,  unconfined 
compression)  stress-strain  curve  in  Figure  2.1b 

Young's  modulus  of  elasticity,  corresponding  to  locked  clas- 
tic region  in  Figure  2.1,  psf  = 3G 

Internal  energy 

In  situ  modulus  of  deformation 
Intact  modulus  of  deformation 

Strain-hardening  modulus,  corresponding  to  locked  plastic 
region  in  Figure  2.1,  psf  = 3G^ 

Tangent  modulus  at  50  percent  ultimate  strength 

Earth  penetrating  weapon 


f Fracture  frequency 

f(c)  Ogive  nose  factor 

g Acceleration  of  gravity 

G Elastic  shear  modulus 


Gtrain-hardening  shear  modulus 

H High  strength 

K Normal  joint  stiffness 

n 
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I Projectile  body  length 

L Low  strength 

m Projectile  mass 

M Medium  strength 

pg  E[l  - exp(  -33)  ] “ f Y in  6 + |f  ^Et  " 7)  Etn 

P Maximum  depth  of  penetration 

P„  Penetra* ion  into  high-strength  rock 

H 

?L  Penetration  into  low-strength  rock 

P„  Penetration  into  medium-strength  rock 

M 

q Material  constant  relating  shock  velocity  to  sound  and 
particle  velocities 

R Radius;  also  Schmidt  hardness 

RQD  Rock  Quality  Designation 

s Fracture  spacing 

t Duration 

V Velocity 

Vf  In  situ  compressional  wave  velocity 

VL  Laboratory  (intact)  compressional  wave  velocity 
Vq  Impact  velocity 

W Weight 

W/A  Sectional  pressure 

Y Yield  strength  of  target  material,  psf 

z Direction  of  penetration 


Y Dry  unit  weight 
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Yq  Gr line i sen  constant 

P 

6 1 - — exp  (-38) 

PP 

At  Time  increment 
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AV, 


Az, 


Incremental  velocity  change  in  the  k^*1  time  step  as  com- 
puted from  Equation  2.1 

Distance  penetrated  during  the  kth  time  step 


3 

2 CFH 


oo 


n=l  ri 


also  1 


P 


♦ (e) 


Density 

Initial  density  of  target  material,  slugs/cu  ft 
Po  exp(ep) 

Uniaxial  compressive  strength 
Limiting  cone  half-angle 
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